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ABSTRACT 

The  multi-loop  Lead  Anfle  servo  loop  is  utilized  as  a  analog  computing 
loop  that  mechanizes  the  range-bearing  torpedo  fire  control  equation. 
This  thesis  applies  several  new  servo  analvsis  techniques,  developes 
system  component  transfer  functions,  analyzes  individual  loops,  and 
examines  qualitative  stability  and  response  of  the  overall  system. 
No  response  or  error  criteria  was  specified.  Based  on  the  author' 3 
experience  and  advice  from  submariners  a  minimum  error  criteria  is 
derived.  This  thesis  preposes  several  modifications  of  the  oripional 
svstem  in  order  to  improve  performance,  points  out  several  deficiencies, 
and  recommends  areas  of  further  study. 
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and  cooperation  extended  to  him  by  the  personnel  he  had  the  pleasure 
of  working  with  at  the  Librascope  Division  of  General  Precision, 
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TABLF  OF  SYMBOLS  AND  ABBREVIATIONS 


Symbol 

B 


Name 


Relative  Target  Bearing 


B-Bg 


Missile  Azimuth 


Bb6 


Torpedo  Lead  Angle 


e(Bb6) 
Bg 

Bms 


Frror  in  Bb6 
Relative  Tube  Bearing 

Torpedo  Angle 


Definition 

When  used  alone,  3  denotes 
the  relative  bearing  to  the 
present  position  of  the  target 
and  is  measured  in  the  Master 
own  srip  horizon  system. 

The  difference  between 
Relative  Target  Bearing  and 
Relative  Tube  Bearirg,  mea- 
sured in  the  horizontal  plane 
of  the  own  ship  horizon  system. 

The  angle  formed  by  the 
intersection  of  the  line  from 
own  ship's  present  position 
to  target's  present  position 
with  the  line  from  torpedo's 
pseudo  or  real  position  to 
target's  future  position. 
Positive  angles  are  measured 
clockwise  from  the  line  connect- 
ing own  ship  and  target,  i.e., 
the  present  relative  bearing 
line.     Prior  to  event  h,  and 
for  all  Preset  problems, 
point  h  is  the  pseudo  position. 
For  Wire  Guided  problems, 
after  event  h,  the  dead  reckoned 
torpedo  position  is  used. 


Relative  bearing  of  the  torpedo 
tube  measured  in  an  own  ship 
horizon  system. 

Angle  between  the  horizontal 
projection  of  the  missile  or 
torpedo  speed  vector  and  the 
present  bearing  line  from  own 
ship  to  the  missile  or  torpedo, 
measured  in  the  horizontal 
plane  of  a  missile  or  torpedo 
horizon  system. 
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S-mbol 


Name 


Definition 


3ts 


Ta-get  Angle 


Co 


Cm 


Own  Ship  Course 
Torpedo  Course 


j(Cm) 

Computed  Correction 
to  Torpedo  Course 

Cmh 

Torpedo  Co  .rse  at  Fvent 
h 

Cmoii 

Relative  Torpedo  Course 
at  Event  h 

ct 


D3y 


DCo 


DCt 


DKho 


Target  Course 


Bearing  Rate 


Own  Ship  Turning  Rate 


Target  Turning  Rate 


Own  Ship  Speed 


Angle  between  the  target  speed 
vector  and  the  present  bearing 
line  from  own  ship  to  target, 
measured  in  the  target  horizon 
system. 

Course  of  own  ship  measured 
in. the  Master  system. 

Course  of  the  torpedo  or 
missile  measured  in  the 
horizontal  plane  of  a  north 
oriented  horizon  system  at 
The  torpedo  or  missile. 


The  course  of  the  torpedo 
measured  in  the  horizontal 
plane  of  a  north  oriented 
horizon  system  at  point  In 

The  course  of  the  torpedo 
measured  in  the  horizontal 
plane  of  an  own  ship  oriented 
horizon  system  at  point  In 

Course  of  the  target  measured 
in  the  target  north  oriented 
horizon  system. 

The  time  rate  of  change  of 
true  target  bearing. 

The  rate  of  change  of  own 
ship  course  with  respect 
to  real  time. 

The  rate  of  change  of  target 
course  with  respect  to  time. 

The  time  rate  of  change  of 
the  distance  traversed  by 
own  ship,  measured  in  the 
horizontal  plane  of  the 
Master  own  ship  horizon 
system. 
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P;  mbol 
PMht 


Name 


Target  Speed 


Fc 


Fo 


G 


•  (G) 

sq(G) 

G6 
Hm 

c(Hm) 
ec(Hnt) 
q(Hm) 
sj(Hm) 


Frequency,  cycles  per 
second 

Closed  loop  transfer 
function 

Open  loop  transfer 
function 

Gyro  Angle 


Error  in  Gyro  Angle 
Angular  Spread 

Gyro  Angle  at  Event  6 
Lenpth  of  Torpedo  Path 


Computed  Length  of 
Torpedo  Path 

Error  in  Computed  Length 
of  Torpedo  Path 

Run  to  Burst  Command 


Correction  to  Torpedo 
Path  Length 


Definition 

The  time  rate  of  change  of 
the  distance  traversed  by 
the  target,  measured  in  the 
horizontal  plane  of  the  target 
horizon  system. 


Angle  betveen  vertical 
plane  thru  line  of  fire  and 
vertical  plane  thru  desired 
missile  centerline  measured 
in  horizontal  plane.  (When 
used  alone  only.) 


A  manual  correction  to  gyro 
angle  to  initiate  an  angular 
spread  of  torpedo  paths. 

The  torpedo  gyro  angle  at 
event  6. 

Actual  distance  traveled  by 
torpedo  from  own  ship  to 
target. 


Manual  correction  to  torpedo 
path  length, 

A  bias  added  to  torpedo  path 
length  to  compensate  for 
torpedo  acceleration  to 
required  speed. 


ix 


Symbol 

sq(Hm) 

Hm26 


Hro6 


Hm7 


Hvm 


Hvmg 


Name 


Enabling  Run  Offset 


Total  Run  to  Burst 


Di  stance  from  Torpedo  to 
Point  6 


Remaining  Run  to  Burst 


Torpedo  Run  Depth 


Depth  Difference 


Definition 

Correction  used  to  shift 
enabling  point  along  the 
torpedo  track. 

The  total  theoretical  length 
of  the  torpedo  path  from  event 
2  to  event  6. 
(Obsolete  symbol  Hmt) 

The  straight  line  distance 
from  the  torpedo's  present 
position  to  point  6. 

The  total  distance  the  torpedo 
must  run  from  its  present  t> 
position  to  its  bursting  position 
(point  7). 

Vertical  distance  the  torpedo 
is  set  to  run  below  the  surface 
of  the  water. 

Vertical  distance  between  the 
torpedo  tube  and  the  missile 
or  torpedo.   (Hvm  -  Hvg) 


K27 

Distance  from  event 
2  to  event  7 

H67 

Linear  Spread  Along  Tar 

get  Track 

J 

Ly2 

Latitude  of  Firing 

Lyo 

Latitude  of  Own  Ship 

Lyp 

Latitude  of  Proofing 

Lyt 

Latitude  of  Target 

k 

One  Thousand  ohms 

K 

Gain 

The  distance  along  the  target  track 
between  the  future  target  position 
(point  6)  and  the  advance  point  (point7) 


The  latitude  of  own  ship  at  the  time 
of  fire. 


The  latitude  at  which  the  torpedo 
was  proofed. 


Symbol 
Mho 

Pdn 


Name 

Own  Ship  Horizontal 
Notion 


Athwart ship  Parallax 
Displacement 


Pdo 


Centerline  Parallax 
Displacement 


Definition 

The  horizontal  distance  made  good 
by  own  ship  during  a  given  inter- 
val of  tine. 

The  component,  in  the  deck  re- 
ference plane  and  perpendicular 
to  own  ship  centerline,  of  the  total 
distance  from  the  origin  of  the 
Master  system  to  the  origin  of  the 
weapon  system.   (The  center  of  the 
torpedo  gyro  or  the  missile  platform. ) 
Pdn  is  measured  in  the  Master  own 
ship  deck  system. 

The  component,  along  own  ship  center- 
line,  of  the  total  distance  from  the 
origin  of  the  Master  system  to  the 
origin  of  the  weapon  system.   (The 
center  of  the  torpedo  gyro.,  or  the 
missile  platform.)  Pdo  is  measured 
in  the  Master  own  ship  deck  system. 


rx,r2< 

R 


R^ots  of  equation 
Present  Slant  Range 


j(R) 
Rh 


Computed  Correction  to 
Range 

Horizontal  Range 


c(Rh) 
j(Rh) 

Rh2 


Computed  Horizontal  R;'nge 
Horizontal  Ranpe  Correction 


Horizontal  Range  to  Target 
at  Event  2 


R,  when  used  alone,  is  the 
present  distance  between 
own  ship  and  target,  mea- 
sured along  a  straight  line 
in  the  Master  system. 


The  present  horizontal  rang© 
from  own  ship  to  the  target's 
present  position  as  seen  by 
radar  or  as  computed  from 
sonar  data,  measured  in  the 
Master  system. 


Computed  correction  to 
horizontal  range. 

The  range,  measured  in 
the  Master  system,  from 
own  ship  to  target  at  time 
of  fire  (event  2). 


Name 


Definition 


Horizontal  R^nge  to    ?each 

Point 


The  range,   measured  in 
the  Mastpr  system,   from 
own  ship  present  position 
to  the  torpedo  reach  point, 
ooint  3. 


t 

Time 

/* 

lT%r 

€ 

Base  of  Natural  Logarithm 

£ 

Frror 

e 

Angular  displacement,  "-radians 

ft; 

System  Input 

Go 

System  Output 

Tf 

Time  of  Flight 

7 

T26 


Urn 


Angular  Velocity,  radians 
per  second 

Time  Constant 

Microfarads 

Time  of  Travel 


Time  to  Event  x  in 
Fig.  1. 

Torpedo  or  Missile  Speed 


j(Um) 


Correction  to  Torpedo 

Speed 


The  total  time  from  time 
of  fire  (event2)  to  time  of 
torpedo  burst  (event  6). 


The  speed  of  the  torpedo  or 
missile  with  respect  to  the 
medium  through  which  it  is 
moving.  The  time  rate  of 
change  of  the  distance  over 
which  it  is  moving. 

A  computed  correction  to 
torpedo  speed. 
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Symbol 
Ym 


V.arre 
Torpedo  Turn  Radius 


Target  Half -Length 


Definition 

The  mean  radius  of  curvature 
of  the  path  taken  by  the  torpedo 
(essentially  circular)  in 
traveling  from  the  reach  point 
(point  3)  to  the  gyro  involute 
point  (point  h) . 

Distance  from  the  stern  of  the 
target  to  the  middle  of  the  target 

(MOT). 
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T.   SUMMARY 

An  analysis  of  individual  loops  was  conducted  and  a  representative 
loop  was  selected  and  presented  herein.  The  individual  loop  as 
designed  will  ooerate  satisfactorily;  however,  using  optimum  settings 
of  input  and  feedback  resistances  determined  in  section  5-1  the  response 
can  be  improved  37.5  per  cent  over  the  best  response  obtained  by  trial 
and  error,  using  settling  time  as  a  criteria.  A  further  improvement 
was  obtained  by  eliminating  the  integral-damping  network;  however, 
considerable  noise  was  incountered  which  would  not  show  up  in  the  output 
but  would  cause  gear  wear  in  the  chain  near  the  servo  motor.  By  chang- 
ing capacitors  C-3  and  C-U  in  figure  lu  to  .25  microfarads  instead 
of  the  presently  installed  1  -nicrofarad,  the  same  response  could  be 
obtained  without  the  undesirable  oerformar.ee  at  steady  state.  The  over- 
all improvement  in  response  using  the  resistances  determined  and  the 
change  in  the  capacitance  was  5b. 2  per  cent. 

The  overall  system  was  considered  but  because  of  the  non-linear  nature 
of  the  resolver  in  the  circuit  exact  response  cculd  not  be  predicted  in 
the  time  allowed  for  analysis.  Linearizing  assumptions  were  made  and  a 
stability  analysis  was  performed.  It  was  found  that  the  system  was  high- 
ly stable  with  a  gain  margin  of  2k   db  and  a  phase  margin  of  72  degrees. 
During  the  analysis  it  was  discovered  that  the  system  was  type  zero 
which  has  an  infinite  steady  state  error  to  a  ramp  input.  Within  the 
area  of  expected  performance  of  the  system  this  error  did  not  appear  to 
be  significant;  however,  if  large  ramp  inputs  are  incountered  or  long 
c  mputing  times  the  error  can  become  significant.  The  system  is  easily 
compensated  to  become  type  one  by  an  additional  amplifier-motor  loop. 


If  it  is  desired  to  use  the  system  as  designed  it  appears  that  by 
careful  scaling  that  the  error  can  be  reduced  significantly  without 
further  compensation. 


CI.    NTiTDUCTION 

2.1  The  Gyro  Lead  Angle  Computing  Loop 

The  gyro  lead  angle  computing  circuit  functions  as  a  part  of  the 
Attack  Director,  Mk  75.  It  is  considered  advantageous  to  famil^a- 
the  reader  with  the  overall  svstem  Drior  to  discussing  the  specific 
circuit  analyzed. 

The  Attack  Director,  Mk  75  performs  as  a  torpedo  data  computer  and  nosit- 
ion  keeper.   Continuous  inputs  of  ranre  and  bearing  information  ai 
from  various  selected  data  rowers. 

As  a  position  keeper  the  Attack  director  renerates  range  and  beai 
tarret  position  for  transmission  to  sonar  via  the  Analyzer  Console, 
from  the  most  recent  values  of  target  course,  soeed,  rar.ge,  and  bear 
or  venerates  range  and  bearing  of  apparent  tarret  position  for  dir< 
ission  to  sonar,  from  the  most  recent  values  cf  automatically 
nanually  inserted  inputs  of  tarret  course  and  soeed. 

As  an  angle  solver  the  Attack  Director  computes  gvro  angle,  enal 

•  lectrically  set  torpedoe  . 

quired  information  ^or  the       1  of  wire  ru  led  torpedoes 
bo  torpedo  course  and  run  to  burst,  and  a  c 
is  aco  -  Lishi  r. 

■  - 
the  bearinr  line  to  a  t     .    rpedo  inputs  ar< 
or  to  time  of  fir&  for  w      .  ided  and  pre 
:sition  is  displave  f  fire  for  t] 


The  Attack  Director,  Mk  75  functions  in  both  automatic  and  iisplay  -"odes. 

The  specific  loop  which  is  the  subject  of  this  analysis  is  a  multi-loop, 
electro-mechanical  servo-mechanism.  The  loop  is  composed  of  eight  individ- 
ual loons,  which  are  interdependent  in  the  operation  of  the  Attack  Director, 
Mk  75.   Listed  below  are  the  elenent  loops  of  the  comoensation  of  lead 
anple  loop  and  their  phvsical  location  in  the  director: 

1.  Gyro  Angle  -  located  in  the  Angle  Solver  Section 

2.  Enabling  Run  -  located  in  the  Angle  Solver  Section 

3.  Torpedo  Lead  Angle  -  located  in  the  Position  Keeper  Display 
Section  * 

li.  Relative  Tarpet  Bearing  -  located  in  the  Position  Keeper 
Display  Section  of  the  Attack  Director 

5.  Relative  Tube  Bearing  -  located  in  the  Ballistic  Computer 
Section 

6.  Weapon  Time  to  Furst  -  located  in  the  Ballistic  Computer 
Section  of  the  Attack  Director 

7.  Target  Speed  -  located  in  the  Position  Keeper  Display 
Section 

8.  Target  Anrle  -  located  in  the  Position  Keeper  Displav 
Section 

The  commutation  of  the  lead  angle  loop  is  a  seauential  computation  loop. 
This  loop  has  four  possible  mechanical  inputs  and  five  possible  electrical 
inputs.  The  outputs  from  this  loop  are  three  dial  readings.   These 
readings  are  valid  onlv  when  the  loop  has  reached  a  steady  state  condit- 
ion. 

In  order  to  explain  the  operation  of  the  computation  of  lead  angle  loop 
several  assumptions  were  made.   Initial  conditions  are  such  that  Gyro 
Angle  |GJ  ,  Torpedo  Lead  Angle  [&b  L] ,   Enabling  Run  [C  (Mm)  j  ,  and 


Weapon  Time  to  Burst  [T*  /,  are  initially  at  some  particular  value.  Another 
condition  is  that  electric  rower  to  all  loops  is  not  applied  until  such 
time  as  the  computation  begins. 

Another  assumption  made  is  that  all  electrical  inputs  are  constants  with 
the  exception  of  the  following  variables. 

1.  Ym'    Sin  Bb6 

2.  Cos  Bts   (Dmht   -Kl    .   T26) 
"*.      Ym'    Cos   Bb6 

)..      Sin  Bts   (Dmht   .    Kl   .    T26) 

With  the  above  in  mind  let  us  consider  the  solution  to  the  problem: 

■» 

The  firing  of  a  preset  torpedo  reouires  solvinp  for  a  predetermined  point 
of  impact  between  torredo  and  the  target,  noted  as  point  6.   The  torpedo 
inout  requirements  are:  run  to  the  enabling  point,  c(Hm);  the  gyro  anrle 
order,  G;  run  dr-pth,  Hvm;  and  a  26  Volt  mark  signal  when  0  <T<180,  for 
automatic  horizontal  search. 

Computations  are  made  to  resolve  into  mutually  perpendicular  components 
the  final  vector  (lp»g)  of  the  torpedo  run  to  the  tarret.  This  is  the 
straight  line  connecting  event  h   with  event  6.   The  components  are  JETY, 
which  is  parallel  to  the  line  of  sight  (LOS)  and  2X,  which  is  perpendicular 
to  it,  (reference  figure  1). 

The  ouantities  2"  X  and  XY  are  predicted  by  two  isolation  amplifiers  summing 
the  applicable  inputs.  These  two  amplifiers  are  the  start  of  a  seouential 
computation.  The  amplifier  inputs  involve  the  parallax  quantities  of  the 
tube,  Pdo  and  Pdn,  the  reach  of  the  torpedo  Rh3,  the  torpedo  turn  radius 
Ym1 ,  the  relative  Target  Bearing  B,  corrected  for  Relative  Tube  Bearing 
Bg,  the  Torpedo  Lead  Angle  Bb6,  the  Target  Angle  Bts,  target  speed  times 
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the  total  run  time  of  the  torpedo   (DMht   .    T?6),    and  a  linear  scread  K67, 
quantity  if  present.      These  quantities  are  related  to  £X  and  HY  bv  the 
following  equations. 

X  =  sin  B    (Pdo+  Rh3   cos  Bp  -  Ym'    sin  Br) 

-  cos  B    (Pdn  +  Rh3   sin  Bp  -  Ym'    cos  Bg) 

+Ym'    cos  Bb6   +   sin  Bts   (OMht    .    K      .    T26   +  H67) 
Y  =  Rh  -   cos   :3(Pdo   +  Rh3   cos  Bp  -  Ym'    cos  Bg) 

-  sin  B    (Pdn  +  Rh3   sin  Bp  -  Ym'    cos  Bg) 

-  Tm'    sin  Bb6  -  cos  Bts   (DMht   .    ^    .   T  26  +  H67) 

The  cuantities^X  and  i-Y  are  inputs  to  the  two  stator  windinps  of  a  resol- 
ver.     A  servo  loop  positions   the  mechanical  inout  of  the  resolver  until 
one  of  the  rotor  outputs  approaches  null.      The  output  of  the  other  rotor 
is  the  vector  sum  of  £X  and  £Y  Hm6.      The  shaft  position  of  the  servo  loop 
is  an  anple  equal  to  the  tan"'  ~£    ■  Bb6,    the  torpedo  lead  ar.ple   (refer- 
ence  fipure  1). 

A  perturbation,  which    includes  a  co-npensation  for  drift  and  gyro  uncag- 
ing time,   is  summed  with  the  null  signal  e(Bb6)   from  the  resolver  in  an 
isolation  a^Dlifier.      The  generation  of  this  perturbation  will  be  explained 
later.     The  output  of  this  isolation  amplifier  is  sent  to  a  variable  gain 
amplifier, VGA,  which  divides  this  sirnal  by   the  other  output  of  the  lead 
angle  resolver. 

The  output  of  this  amplifier  is  the  error  in  Bb6. 
C(B6<J  =  jiZtCos  SbC-TV  S,+0£6j   /•  ^~^  (S/as  ly/0 

-  $"  //a  J  -  Oc0  ry]  c(//*,jj    -r  /V/»  c 
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This  error  is  computed  by  the  relation  9  -^where  0  is  the  error  in  Bb6, 
s  is  an  error  in  the  difference  in  length  of  a  particular  component  of 
£  X  and  £  Y  and  the  perturbation  mentioned  above,  and  r  is  the  vector 
sum  of  IX  and  Z.  Y,  Hm6. 

The  vector  sum  of  IE  X  and  £T  is  the  length  of  torpedo  travel  from  point 
k   to  point  6.  This  signal  is  sent  another  isolation  amplifier,  whose 
function  will  be  explained  later.  The  3haft  position  of  torpedo  lead 
angle  (Bb6)  is  displayed  on  the  outer  dial  of  the  "Target  Group"  and  "Own 
Ship  Group"  dial  display  in  the  Position  Keeper  Display  Section. 

The  Bb6  shaft  position  is  added  to  he  target  angle  (Bts)  by  a  mechanical 
differential  (Bb6  +  Bts  ■  I).  The  output  of  this  differential  positions 
a  limit  cam,  which  operates  a  micro  switch.  The  normally  closed  position 
of  the  microswitch  gives  a  right  turn  mark  signal  to  the  torpedo  for 
automatic  horizontal  search,  NC,  0<I<.180  •  The  normally  open  position 
provides  no  mark  signal  and  therefore,  aleft  turn  is  in  order.  This  turn 
signal  is  available  as  a  oreset  torpedo  function  that  is  activated  at 
the  enabling  point. 

The  Bb6  loop  also  positions  a  two  speed  control  transmitter  synchro  sy- 
stem.  These  CX's  send  an  electrical  analog  signal  proportional  to  the 
Torpedo  Lead  Angle,  Bb6,  to  a  two  apeed  control  differential  transmitter 
svnchro  system  that  is  mechaically  referenced  with  the  Relative  Target 
Bearing,  B.  These  CDX's  add  the  two  functions  and  yield  Relative  Torpedo 
Course,  Cmoh.  The  CDX's  then  send  this  signal  to  the  Own  Ship's  course 
modulo.  The  Cmoli  signal  al30  references  a  two  apeed  control  differential 
transmitter  synchro  system  that  is  mechanically  referenced  with  Relative 
Tube  Bearing,  Bg.  The  output  of  the  relative  tube  bearing  CDX's  reference 
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a  two  speed  CDX  svstem  on  the  Angular  Spread  module  sq  (G).  This  ^odule 
is  located  on  the  firing  panel.   The  output  of  these  CDX' s  reference  a 
two  soeed  control  transformer  svstem  on  the  Gyro  Angle  Order  module,  G. 
These  CT's  are  the  feedback  element  for  the  loop.  The  eauation  involved 
here  is:  [q  =  Cr*o4-  Bu      +^(<?j]  . 

Sach  of  the  terms  on  the  r^  rht  side  of  the  above  eouation  for  G  have  been 
added  by  the  use  of  two  soeed  svnchro  systems  (CDy).  This  signal  refer- 
ences the  gyro  angle  module  feedback  CT's.  These  CT's  have  a  shaft  position 
proportional  to  Gyro  Angle,  G,  and  an  output,  which  through  a  servo  amp- 
lifier and  motor,  positions  the  pvro  angle  order  shaft. 

The  gyro  angle  order  gear  train  positions  the  wioer  on  a  DC  potentiometer 
(rot).  This  G  pot  is  referenced  with  TJTtimes  torpedo  right  t'.;rn  radius, 
TT"  VfH      anc^  ~^~  times  torpedo  left  turn  radus,  77""  Vm*      .  The  center 
tap  of  the  pot  is  grounded.   The  voltage  from  the  not  is  positive  for  both 
ripht  and  left  turns,  but  it  is  not  of  the  sa^e  magnitude.  The  voltage 
output  at  the  wiper  position  of  this  pot  is  then  proportional  to  the  wiper 
position,  a  function  of  G,  and  the  reference  voltare.   The  analog  signal 

is  then  the  product  of  these  two  inruts;  hence  the  arc  distance  traveled 

77- 

about  the  gvro  angle,  Ym'  .  G  .— -,.  This  signal  is  sent  to  an  isolation 

amplifier,  figure  2.  Other  inputs  to  this  amplifier  are:  The  distance 

to  the  reach  point,  Rh3,  from  the  ballistic  constants  network;  the  vector 

sum,  £X  X  sin  Bb6  +  X   Y  cos  Bb6  -  Hm6J ,  from  the  lead  angle  resolver 

previously  discussed;  and  a  correction  for  torpedo  run  depth  Hvm,  and 

own  ship  depth  Hvo,  Hvmg   ,  from  the  ballistic  constants  network. 

Kn3 

The  output  of  this  amplifier  is  proportional  to  Hm26.  This  voltage  is 
sent  to  an  impedance  matching  transformer  in  the  input  network  of  a  servo 
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amplifier.  The  output  of  this  network  is  the  output  of  the  feedback 
vernistat  Hm26  on  the  Enabling  Run  servo  module.  The  amplifier  drives 
the  servo  motor  until  the  two  inputs  are  equal  in  magnitude  and  180  out 
of  phase.  The  motor  has  then  driven  out  a  mechanical  analog  position 
proportional  to  Hm26.  This  shaft  position  may  be  added  to  or  subtracted 
from  manually  by  a  hand  crank  input  labeled  Enabling  Run  Offset,  sq(Hm). 
This  is  accomplished  by  the  use  of  a  mechanical  differential.  The  amount 
of  offset  is  recorded  on  a  counter  for  display. 

The  Enabling  Run  module  visually  displays  computed  Enabling  Run,  c(Hm), 
which  is  the  summation  of  Hm26  and^the  offset  sq(Hm).  One  vernistat  on 
this  loop  sends  a  voltage  proportional  to  the  computed  Enabling  Run  c(Hm), 
to  the  Run  to  Burst  module  as  a  reference  source  for  a  vernistat  feedback 
loop.  This  particular  total  run  to  burst,  (_Hm26  w/g  -  NR  ■  c(Hm)  w/g  + 
q(Hm)J  ,  is  used  only  in  the  wire  guided  solution  and  will  not  be  discussed 
here. 

A  separate  cam  for  each  weapon  on  the  computed  Enabling  Run  shaft,  c(Hm) 
figure  3»  sends  a  within  limit  signal  to  the  gyro  angle  limit  cam.  If 
the  gyro  limit  cams  are  within  limits,  then  they  send  a  signal  to  the 
depth  stratum  limit  cams.  The  run  depth  cam,  then  sends  a  signal  to  the 
within  limit  circuitry  in  the  weapon  monitor  unit. 

The  Enabling  Run  module  position  six  one  speed  control  transformers  that 
send  an  analog  error  signal,  ec(Hm),  to  a  servo  motor  in  the  torpedo. 
There  is  a  control  transmitter,  mechanically  linked  to  a  servo  motor  in 
each  torpedo.  This  CX  feeds  a  follow-up  signal  oc(Hm)  back  to  the  CT's 
in  the  Enabling  Run  module  to  comolete  the  servo  loop. 

The  second  vernistat  in  the  Enabling  Run  module  is  referenced  with  the 
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difference  between  the  sine  of  firing  latitude  and  the  sine  of  proof  i  ng 
latitude  times  half  the  product  of  the  enabling  run  time  and  target  turning 
rate,  this  m.  nus  the  signal  for  gyro  uncaging  time,  times  own  ship  turn- 
inp  rate.  These  two  quantities  are  multiplied  by  enablinp  run  by  the  use 
of  the  vernistat. 

The  output  of  this  vornistat  is  sent  as  the  perturbation  to  the  lead  anrle 
servo  loop,  figure  U,  previously  discussed. 

The  analog  voltage,  c(Hm),  from  the  Enabling  Run  Vernistat,  is  sent  to 
the  total  run  to  hurst  module,  figure  5,  and  to  an  isolation  apip'ifier, 
which  sums  computed  Enabling  Run,  c(Hm),  and  a  torpedo  oath  length  correct- 
ion, Hmcorr 

Hm     »  -  Sg(Hm)  ±  Hvmp  (Urn)   +  Hvmfr   (Hm) 
corr.     L  c  J 

(Note:  The  Hvmr  /,,  x  term  is  used  only  with  the  Mk  37  Mod  0,  high  speed 
torpedoes.)  The  output  of  this  amplifier  is  total  torpedo  run,  (-)  Km27. 
This  signal  is  sent  to  an  impedance  matching  transformer  in  the  in- ut 
network  of  the  isolation  amplifier  of  the  Weapon  Time  to  Burst  servo  system. 
A  vornistat  on  the  Weapon  Time  to  Burst  module  is  the  other  inrut  to  the 
amplifier.  This  vernistat  is  referenced  with  toroedo  or  missile  speed,  Urn', 
from  an  isolation  amplifier  having  3  ballistic  constants  f(Hmv),  Ut,  i(Um) 
for  inouts.  The  f(Hvm)  signal,  a  function  of  depth,  is  a  correction  to  the 
torpedo  or  missile  speed,  Ura.  This  compensation  is  generated  on  the  Run 
Depth  module  by  a  pot  positioned  proportionally  to  Run  Depth,  Hvm,  and 
referenced  from  the  ballistic  constants  network  with  a  function  f(Hvm). 
The  voltage  output  of  the  T26  computing  vernistat  is  proportional  to 
torpedo  speed  Urn1  times  weapon  time  to  burst  T26.  When  the  product  of 
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of  speed  and  time  (ecual  to  total  run  distance  Hm-t),  is  compared  with 
(  -  )H?7  and  passed  throuph  an  isolation  amplifier,  it  is  an  error  sipnal, 
e(T26),  used  to  position  the  Weapon  Time  to  Burst,  T26,  servo  system 
which  displavs  T26.   A  vernistat  referenced  with  the  difference  between 
the  sine  of  proof inp  latitude  and  the  sine  of  firinp  latitude  time3  half 
the  torpodo  turning  rate  and  positioned  by  the  T26  servo  system  has  an 

analop  voltape  output  equal  to  the  Toduct  of  v.eapon  tine  to  burst  and  the 

fDCm      .    .  .1 

condensation  for  firinp  latitude  I  ~T~  •  T26  \sir\   Lyp  -  sin  Lv2)J  .   This 

sipnal  is  then  summed  at  the  input  of  an  isolation  amplifier  with  a  con- 
tribution from  the  ballistic  constant  network.   ( -DCo  .  Td),  where  DCo 
is  the  Own  Ship  Turning  Rate  and  Td  is  the  Gvro  Uncapinp  Time.  The  out- 
out  of  the  isolation  amplifier  is  used  to  reference  the  vernistat  on  the 
Enablinp  Run  nodule  previously  discussed. 

Another  vernistat  in  the  Weapon  Time  to  Burst  module  sends  a  sipnal 
prooortional  to  T26  to  reference  a  vernistat  on  the  Tarpet  Speed  module. 

This  vernistat' s  output  is  a  multiple  of  Tarpet  Speed,  DMht;  Weapon  Time 
to  Burst,  T26  and  K^,  to  vield  the  distance  the  tarpet  travels  from  time 
of  fire  to  time  of  impact.  This  sipnal  is  sent  to  an  isolation  amplifier, 
figure  2,  that  sums  this  distance  and  linear  spread  quantity  H67,  if 
present.   The  linear  spread  inout  is  located  in  the  firinp  nanel.  Linear 
spread  is  never  used  concurrently  with  Angular  Spread  auantities,  sq(G), 
previously  discussed.  The  output  of  this  amplifier,  the  total  distance 
travelled,  is  [(Tf  .  DMht  .  Kx)  ♦  H67]   or  [(T26  .  DMht  .  Kx)  ♦  H67J  , 
and  references  a  resolver  in  the  Tarpet  Anple  servo  system.  The  output 
of  this  resolver  is  the  sine  and  cosine  of  Tarret  Anple,  Bts,  times  the 
reference  signal.  These  two  signals  become  part  of  the  inputs  to  the 
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two  isolation  amplifiers  which  vield  5"  X  and  IT, 

Generation  of  one  of  the  other  inputs  to  each  of  these  amplifiers  bepins 
with  the  output  of  two  isolation  amplifiers.  The  input  to  these  amplifi- 
ers is  dependent  uoon  the  value  of  the  Gyro  Angle,  C.  Three  conditions 
can  be  considered: 

(1)  For  170°^  G  *10?  and  3$0°>  G  >190*;  the  inputs  into  the  turn- 
radi us-isolation  amplifier  is  just  Tm'  and  the  incut  into  the 
reach-point-isolation  amplifier  is  Rh3  ♦  Hvmg™~  where  Hvmgpj- 
is  the  correction  to  Rh3  for  depth  difference,  (Hvo  -  Ph  -  Hvm). 
Under  these  conditions  X-OC  and  £Y  are  comouted  as: 

II  «  [Pdo  sin  B  -  Pdn  cos  B  ♦  (Rh3  ♦  Hvr^Rh3^  sin  (B  "  bp)  "vm' 
cos  (B  -  Bp)  +  Ym'  cos  Bb6  +  (DMht  .  T26  .  Kx  +  H67)  sin  BtsJ 

(2)  For  350<G<10a  the  inrut  to  the  pvro-turr-radius-isolation 
amplifier  is  zero;  however,  the  reach-noint-isolation  amplifier 

's  r.ow  Ph3  ♦  hvmgp^-)  +  (Ym'  .  G  . —   )   .   Tn  addition 

to  this  chanpe,  additional  circuitry  removes  (-Ym"  sin  Bb  . ) 
and  (Ym'  cos  Bb6)  from  the  JE  X  e-nd  JY  isolation  amplifiers. 

Hence  £*  X  and  X  Y  are  commuted  as: 

IX=   Pdo  sin  B  -  Pdn  cos  B  ♦  (Rh3  +  Ym1  .  G  .  JE-     ) 


T   Y  =  Rh  -  Pdo  cos  B  +  Pdn  sin  B  +  (Rh3  ♦  Hvmpnv -,  +  Yir,1  G.  —  ) 
cos  (B  -  Bp)  +  (DMht  .  T26  .  K±   +  H67)  cos  Dts 

(3)  For  l''C  <  G  <  190°  the  innut  to  the  pvro- turn-radius-isolation 
amplifier  is  2  Ym1  or  twice  the  valne  used  in  case  1  above.  The 
reach-ooint-isolation  amplifier  is  the  same  as  in  case  1  and  as 
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before,  additional  circuitry  removes  (Ym'  sin  Bb6)  and  Ym'  cos 
Bb6)  from  the  SIX  and  51 Y  isolation  amplifiers. 

Furthermore,  it  should  be  noted  that  the  quantity  Ym'  used  in 
(l)  is  actually  Ym^  or  Ynv,*,  depending  on  whether  the  condition 
360 >G  >180  exists  or  not. 

In  case  (2)  the  same  thing  is  true,  however,  its  determination 
is  due  to  the  position  of  a  not  positioned  in  the  Gyro  Angle 

module,  figure  6. 

In  case  (3)  Ym'  implies  2Ym'^  or  2Ym'r,  depending  on  how  the  Gyro 
Anple  develops. 

These  sipnals  described  above  for  case  one,  two,  and  three,  reference 
the  stators  on  a  resolver  on  the  Relative  Tube  Bearing  nodule.  This 
module  has  two  positions.  One  position  is  for  the  rearinp  of  the 
port  tube  and  other  position  is  for  the  bearing  of  the  starboard  tube. 
The  input  signals  come  from  the  own  ship  section  where  the.-e  two 
bearings  are  stabilized.   The  distance  to  the  Peach  ;  oint  (point  3) 
and  gyro-turn  radius  are  measured  along  and  associated  with  the  tube 
bearng  module  references  these  ouantities  to  own  ship's  centerline, 
(Own  Ship's  Course,  Co.). 

To  these  quantities  are  added  the  appropriate  parallax  components, 
Pdo,  Pdn,  of  the  ship  by  summing  the  applicable  input  in  two  isolation 
anplifiers.  The  output  of  these  two  amplifiers  then  reference  a  re- 
solver  on  the  Relative  Target  Bearing  module.  This  resolver  rotates 
the  input  quantities  to  distances  along  and  across  the  bearing 
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line  (US).      These  quantities  are  then  sent  to  the  summing  amplifiers 
mentioned  before  and  veld  Z  X  and  £  Y. 

2.2  Cbiectives  of  Analysis 

The  analvsis  was  undertaken  to  determine  the  response  characteristics 
of  the  individual  loops,  the  overall  loop,  check  the  system  for  stab- 
ility, recommend  potentiometer  settings  for  optimum  response  and  if 
found  necessary  recommend  compensation  networks. 

2.3  Period  of  Time  Allowed  for  Analysis 

Work  on  this  thesis  was  accomplished  during  eirht  weeks  of  the  summer 
period  between  the  second  and  third  year  of  the  Ordance  Engineering 
Guided  Missiles  Curriculum  ending  in  December  of  I960. 
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Ill,  Review  or  Introduction  to  Some  Recent  Developements  In  Servo 
Analysis  Techniques 

3.1  Determinantal  Method  of  Analysis.. 

The  theory  of  simple  sinple  loop  control  systems  is  well  developed,  and 
the  block  diagrams  of  such  systems  are  not  complicated.  On  the  other 
hand,  the  block  diagrams  and  theory  of  multi-loop,  multi-input,  multi- 
output,  and  multi-coupled  control  systems  is  not  very  well  developed 
because  of  the  complexity  of  the  problem  involved.  The  block  diagrams 
of  such  systems  are  complicated  and  difficult  to  intepret.  One  approach 
to  the  solution  of  such  a  problem  is  to  formulate  and  solve  a  set  of 
simultanious  equations  for  the  system.  This  can  best  be  accomplished 
by  arranging  such  systems  into  a  standard  block  diagram.  Arranging  the 
system  in  this  manner  greatly  facilitates  the  writing  of  the  equations. 

Any  system  which  has  summing  points  in  the  main  transmission  path  can 
be  arranged  in  a  block  diagram  by  suitable  block  diagram  manipulation 
as  shown  in  figure  7.  Writing  the  equations  at  the  various  nodes: 
At  node  a 

a  "  Ra  -  <Gab-  GAa>a  "  ^W  GBa)b  "  (Gcd-  GC^ 
A  more  convenient  notation  is: 

a  -  Ra  "  <Gaa)a  "  (°ba)b  "  <Gca)c 

At  node 

b  -  Rb  -  (Gbb>b  *  <Gab>a  "  <Gcd)c 

CHU,  T. ,   A  Generalized  Theory  of  Linear  Multi-Loop  Automatic 
Control  Systems,   Doctoral  Dissertation,   MIT,   1953 
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at  node  d 


Rearranging  to  bring  out  the  symmetry 

C-  Oo-kt   d  +(j+v»*^+  (Gc»Jc  +  (G<tbJd   -   ^b 

C"  G* c)  fit      -     {fff)k*d*6ma)C*(G4m)d   '£• 

Now,   solving  by  Cramer's  Rule 


4    - 


(J>^J  Gfca.  '     Geo.  Gda 


-  Ga.4 


-Gb4 


-Gc4  \lfl344j  /**/*>    V/A\,e*J*4. 


For  simplicity  let  us  assume  the  inputs 
Then  the  output  will  be       0  -     (£dc  j   d 


/CJ*/ 


a 
-  Ga.k> 


C 
Gca 


J->G<i*j  Gci 


d 

0 


C  Oj~  *"+!*/ 


D  -  Gde     *     .-  G 


a.d 
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*  /r    tf<*  -  O         J*/0      fib   ^  O 


D  - 


It  is  noted  that  A    is  equal  to  the  characteristic  equation  of  the  sys- 
tem. For  an  example,  take  the  four  node  system  in  figure  8. 
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&    *      /  +    Gz  G3  G*  <»V  +<>*  (  ?*  G*Gi  )    -  characteristic 
equation  of  the  system. 
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0  =     6U*; 


'■»«'»/»/! 


y^&J  /     r    A/  Wi  w>   v^^    r  vjl    ^*   ^*  "->" 

Now  putting  the  transfer  functions  in  the  form  of  numerator  divided  by 
denominator. 


/  +       ki  AS*  A/.  /s+  +~-       j       A/i  //.     '■*  /y~~ 

/      ' r~ 


Which  when  simplified  becomes: 

^     IT 

Ai  ^  ^y  -^j-  *  ^  z^*^*  Vy  ^^  /-  /*<*  ^y  ^y  ^.r 


/r  =  *L  ^k   ^  A/i  A/¥ 


The  roots  of  the  transfer  function  are  the  zeros  of  the  characteristic 
equation. 

Dz  &j    &¥  #s-     *  *'   #*  #*  #*  &*~    +     /**    Ms   AS*   As r    =  O 

Now  if  we  desire  to  determine  the  effect  of  the  feedback  path  G5,  we  can 
rearrange  the  characteristic  equation  in  the  form 

The  quantity  in  the  brackets  can  be  equated  to  zero  and  by  root  locus 
techniques  the  roots  can  be  determined.  The  characteristic  equation  will 
then  be  fis   [    (**Aj(S  t  lj  ^+  Aj^*  **/j   +  N*   N*  N*  N* 


-  0 
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The  final  root  locus  equation  will  then  be: 
The  transfer  function  will  then  be: 


(J* P.J  (S+f*/     ■ 


3 #2  Determination  of  the  Composite  Open  Loop  Transfer  Function  and 
^Telocity  Constant  for  Multi-Loop  Systems 


AS* 

/T  c  />-/;/  =■  p~* 

/*       A    PI   P"*      =   P^ 

c     *.      * 

/  A  /^ 

In  terms  of  numerator  and  denominator  of  the  closed  loop  transfer  function 

f-  ,       *± 

to 

A.  -  /^ 

Now  applying  the  classical  definition  of  the  velocity  constant 

S  +    O 


Jr.  = 


-  -^.  s    ^ 


5^-o    <?c  -  /V 


«. 


It  is  convenient  to  determine  the  open  loop  transfer  function  in  this 
manner  and  then  apply  root  locus  techniques. 
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After  the  open  loop  function  has  been  determined,  it  is  auiet  easy  to 
determine  the  type  and  order  of  the  system.  It  will  have  the  form: 


The  denominator  of  the  open  loop  function  is  unfactored;  the  S  can  be 
factored  out  prior  to  performing  the  root  locus  or  K  versus  real  roots 
thus  the  type  and  number  of  the  system  can  be  determined.  If  the  desired 
type  is  not  obtained,  this  is  a  criteria  for  selection  of  suitable  comp- 
ensation. In  any  event  it  avoids  the  necessity  of  finding  the  roots  of 

the  system  to  determine  the  type  number. 

-. 
It  is  often  necessary  to  perform  a  root  locus  or  K  versus  real  root  plot 

to  determine  the  poles  of  the  open  loop  transfer  function.  As  an  example, 

consider  the  example  of  section  3»1»   The  closed  loop  transfer  function 


is    F  « 


F** dl 

4c  -  >n 


j-  A-,   &*•  "<  ***  » 

/~a      — 


£±    £<3  £¥  £j-   +  SY<l/Ys  /V+.  A^~ 


The  denominator  then  can  be  factored. 

D*  03  Dm  #*   +  ^  ^  **  -*»     ~  a 

/Kt   /v*  /*V  fs    _    (T  J7r 


_  =  e 


^*  03  &M  Os- 
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The  roots  are  then  determined  and  the  open  loop  transfer  then  becomes 


fl   - 


A.    tJr  An  A?*  A/± 


The  closed  loop  roots  are  then  determined  in  the  usual  manner. 

3.3  Root  Locus  Gain  versus  Real  Roots? 

The  usual  method  of  determining  the  roots  of  a  closed  loop  function  is  to 
apply  the  familiar  root  locus  principals.  A  somewhat  similiar  technique 
which  embodies  some  of  the  same  principals  is  the  method  of  plotting 
K  (root  locus  gain)  versus  the  real  roots. 


For  illustration,  let  us  consider  this  unity  feedback  loop  shown  below. 
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The  closed  loop  transfer  function  is 


/£--- 


F. 


O  Ai 


K  = 


The  roots  of  Fc  are  where    0»  t-  A/*  s  o  or 

The  open  loop  transfer  function  has  this  f ormi 


QM) 


04    +  A/* 

0* 


*6J 
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2.  Wheeler,  R.  C.  H.,  Unpublished  Class  Notes 
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Now  plotting  the  poles  and  zeros  on  the  S  -  plane: 


(j~> 


Fig.  9  COMPUTATION  OF  ROOT  LOCUS 
GAIN  FROM  S  PLANE 

By  the  Evan's  method  the  closed  loop  roots  are  determined  by  summing  angles 

to  180°  and  plotting  the  locus  of  the  roots.  Then  knowing  the  gain  associated 

with  the  system,  a  cut  and  try  process  is  performed  to  locate  the  roots  /by  the 


relationship  iC-  L1X±  .  This  process  is  necessary  along  the  real  axis  where 
root  locus  exists.  If  these  values  along  the  axis  are  plotted  the  following 
curves  will  result.  All  the  curves  shown  do  not  result  from  every  open  loop 
transfer  function  but  are  shown  as  an  example  of  all  possible  configeration. 
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Fig.  10  PLOT  OF  ROOT  LOCUS  GAIN 
VERSUS  REAL  ROOTS 

It  can  be  shown  that  the  maximum  between  P.  and  r,   is  the  "oreakaway  point" 

and  the  minimum  between  if,  and  -iAis  the  entrance  point  in  the  Evan*s  Root 

Locus  Method. 

If  the  system  gain  is  given  or  can  be  chosen,  a  horizontal  line  drawn  which 
intersects  these  curves  yields  the  real  roots  of  the  closed  loop  function. 
If  it  does  not  intersect  as  between  /•  and  /^  ,  complex  roots  exist.  Also 
between  zf3  and  /i  ,  one  real  root  and  two  complex  roots  exist  with  the 
gain  chosen  as  shown. 
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This  method  works  equally  well  with  complex  poles  or  zeros  in  the  open  loop 
function  by  using  the  same  relationship  for  gain  /C  —  ~^~  >   hut  one 
must  remember  that  the  distances  to  these  complex  quantities  to  the  point 
at  which  the  gain  is  being  computed  must  be  included. 

From  the  foregoing  discussion,  it  is  seen  that  all  the  real  roots  of  the 
closed  loop  can  be  determined  without  plotting  the  root  locus. 

It  is  a  fact  that  the  denominator  of  the  closed  loop  function  which  is  the 

characteristic  equation  of  the  system  has  the  following  form  and  possesses 

-- 
the  below  listed  qualities: 

Ax*    t  BX*~'  +ck  *" U   =  o 

A  -  1 

B  •  Sum  of  the  roots  of  the  closed  loop  function 

C  •  Product  of  the  roots  of  the  closed  loop  function 

Knowing  the  real  roots  from  the  plot  of  "K  versus  Real  Roots"  and  usinr  the 

relationship   X.  -r  ft*  +   <•  j  *"  fc*  +j  ^c)    *  <P<- ~  j  L**c)    =  & 

^,  +  /i*+A3 +  a.<<  *   y 

we  can  now  solve  for  the  real  part  of  the  complex  pair.  Knowing  this  we  can 
use  the  second  relationship 

x,  fix     ti  (f»c   f  juijC^-y^*>  3  u 

which  becomes 

x.   fix   ij  (  +  *  +  v^J-J  =  U 

Since       itj    w  ±   1 ,  /a/*/  ft  are  known,  we  can  solve  for  Wc  thus  all  roots  are 

known. 
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This  method  provides  a  criteria  for  selection  of  system  gain  since  the  effect 
of  root  location  can  be  readily  seen,  and  is  faster  than  the  rather 
laborious  angle  summation  process  necessary  in  the  Evan's  Method  since  in 
most  cases  only  three  points  must  be  computed  to  determine  a  given  curve. 
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IV.   Derivation  of  Transfer  Functions  of  Components  and  Minim  ;  n 
Velocity  Constant  P'^ouirement 


h.l  Servo  Atpl^  fier 

The  frequency  resoonse  of  the  amplifier  was  measured  and  found  to  be  flat 
to  +  1  db  to  7, COO  kc.   Since  this  is  beyond  the  system  bandwidth  the 
simplification  of  constant  pain  will  be  used  throughout  this  thesis. 
It  was  found  that  this  pain  was  approximately  100,000. 

Because  the  input  sipnal  and  the  various  feedbacks  are  not  summed  at  a 
common  node  (see  Fip.  11),  the  pain  could  not  be  considered  in  the 
conventional  manner.  Rather  the  pain  seen  by  the  various  sirnals  is 
different.   It  was  found  more  convenient  to  consider  the  pain  as  a 
function  of  the  variable  resistors.  This  can  best  be  seen  oy  considerinp 
the  equivalent  circuit  of  Fip.  12. 
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Fip.  12  EQUIVALENT  CIRCUIT  OF  INPUT  NETWORK  TO  SERVO  AMPLIFIER 
The  pains  as  a  function  of  the  resistors  R2,  Rz  and  Ry  is  shown  in  5.1, 
Analysis  of  the  Individual  Loop.  The  entire  closed  loop  eauivalent 
circuit  is  shown  in  fipure  18  and  simplified  block  diapram  equivalent 
is  shown  in  fipure  21. 
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h.2   Tn'^ral  Network 

The  intepral-dampinp  network  is  a  means  of  compensatinp  the  amplifier 
and  is  shown  schematically  in  Fig.  13.  The  equivalent  circuit  is  shown 
in  Fip.  Ik. 


Fig.  1U  EQUIVALENT  CIRCUIT  OF  INTEGRAL-DAW  NETWORK 
From  Fip.  lh  it  can  be  shown  that  the  transfer  function  is: 

Q<>    _  5  i*.  a  j 


Siraol '  fvinp  : 


tf.C.  A>*   C4  a,  J.  Aid 


From  Fip.  13  the  various  values  of  resistance  and  capacitance  were 


obtained  and  substituted  into  the  above  transfer  function  which  veilds 
when  simplified:  "~r*;        3  *  /•  /J-f*      t*//so 

Factoring:  -■  —  —   _  _, 

A  freauency  response  analysis  shown  in  Fip.  15  revealed  that  the  actual 
transfer  function  was  -=2.  -       Z4-C&  s     .  This  transfer  function 
was  used  in  the  analysis. 

It  was  found  later  in  the  analysis  that  a  value  of  Co  ■  C  ■  .2S^/f 
produced  an  improved  response  of  the  individual  loop.  Proceeding  as 

above   the  r^sultinp  transfer  function  is  ___flf_ ►.  It  is  included 

in  this  section  for  convenience. 
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It- 3  Servo  ;totor 

The  Servo  motors  Mk  12,  Mod  0  used  were  manufactured  by  KEARFOTT  COMPANY, 
INC.,  LITTLE  FALLS,  NEW  JERSEY.  To  obtain  the  dynamic  constants,  it  was 
recommended  that  the  following  relationships  be  used. 


/*" 


4  c  *?**>/<**'*     f£~*     *£<. 


Using  a  plot  of  number  of  gear  teeth  versus  inertia  in  OjS-/^*   the 
total  reflected  inertia  was  found  and  added  to  the  given  inertia  of  the 
motor.  The  total  inertia  was  found  to  be  <TS   •>■''<*  X/s>~      J£-/*st'  o* 


J '  *     S-7  8      Gfw-    Cl»i 

/.  o  -     JO 


O  2   -  /A/ 


/^  = 


+  777Z  +  <,  *  *  *700  *"» 


By  suitable  change  of  units  this  becomes: 


/=-  »    z  fo      0y~*  ■  c*t  -   i^*- 


/e/9i> 
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h    =>     .023/    5 /*c. 


Go     a  3T.  5*6 

^  $<>  .  Oi  3/  s    *  /} 

Go  /S« 

When  this  transfer  function  was  used  in  the  computation  to  determine  the 
closed  loop  response  and  checked  against  measured  response,  the  results 
were  found  to  be  35  per  cent  in  error.  Several  different  pole- zero  con- 
figerations  obtained  by  using  various  values  of  tach. meter  feedback  were 
used  and  the  35  per  cent  error  was  found  to  be  consistant.  It  became 
obvious  that  the  motor  time  constant  was  in  error.  With  this  in  view,  an 
investigation  to  determine  the  response  or  a  second  order  system  was  attempt- 
ed by  removing  the  tachometer  and  integral  networks.   Transient  response 
curves  were  obtained  at  various  amplifier  gain  settinps  and  found  to  be 
almost  identical  for  any  setting  of  gain.   The  slight  change  was  accredited 
to  the  pole  of  the  a~plifier  at  60C0.  From  these  curves  it  can  be  seen 
that  ts  (settling  time)  *  .355  U3ing  the  relationship t-~  it  was  fo  md  that 
o^  -    A*./   or  2&  ~  2.Q»2-      .  This  can  be  seen  readily  bv  looking  at 
the  root  locus  plot  of  the  second  order  system  shown  in  figure  16. 
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Fig.  16  ROOT  LOCUS  OF  INDIVIDUAL  LOOP  AS 
A  SECOND  ORDER  SYSTEM 


From  the  response  oC  s       V-  / 

The  motor  transfer  function  then  becomes: 

,  Sao 


3(   S+2S-2J 


As  a  further  check  on  the  accuracy,  a  frequency  response  was  performed.  The 
plot  of  the  magnitude  M  in  decibels  is  shown  in  Fig.  17.  It  was  originally 
intended  to  plot  the  results  on  a  Nichols  chart,  transfer  these  results  to 
a  Bodie  plot  to  find  the  break  point  or  time  constant.  This  was  found  too 
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impractical  because  there  was  enough  backlash  present  to  cause  the  phase 
ar.rle  to  shift  irrationally.  Using  the  value  of  Mpv  from  Fir.  17,  second 
order  curves  were  utilized  to  determine  £     the  dampinr  factor.  Knowing 
£  and  the  natural  frequency  at  which  it  occured,  it  was  found  that 

cL  -  i"  ^  h 

Averaging  the  results  obtained  from  the  transient  and  frequency  response, 

it  was  found  that  the  motor  transfer  function  was: 

-. 

"  Si  S+Z8) 

This  transfer  function  was  used  in  all  the  computations  and  was  found  to 
yeild  exccllant  results. 

[,.3  Derivation  of  Theoretical  'iridium  Kv  Requirement 
The  velocity  constant  reouirement  was  not  specified.  Since  the  system 
under  study  is  a  computation  mechanism  it  was  imoeritive  that  a  realistic 
minimum  value  of  accuracy  required  be  determined  and  used  as  a  roal  in 
the  analysis.  Within  reason  it  would,  of  course,  be  desireable  to  have 
Kv  as  large  as  possible  in  order  to  reduce  steady  state  errors,  but  at 
the  same  time  the  system  need  not  be  designed  to  have  greater  accuracy 
than  sensing  devices  within  the  fire  control  equipment. 

4 

In  the  submarine  fire  control  oroblem,  bearing  rate  in  the  early  stages 
of  the  approach  will  generally  ranre'  between  0-5°  per  minute  (rarelv 
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exceeding  2  per  minute).  As  the  submarine  and  tarpet  closes,  the  bearing 
rate  will  increase,  depending  upon  the  relative  motion  of  the  two  and  their 
ranpe.  At  verv  close  ranres  (1000  yards  or  less)  bearinp  rates  may  become 
as  hiph  as  60  per  minute,  although  this  is  rare  in  typical  problems  en- 
countered by  present-day  hiph  sneed  submarines. 


Consider  the  extreme  case  which  could  conceivably  be  encountered  by  two 
very-high  speed  submarines,  each  travelling  50  knots  on  opposite  courses, 
and  passinp  within  100C  yards  of  one  another.  Their  relative  speed  is 
100  knots  (nautical  miles  per  hour)  and  their  ranpe  is  0.5>  nautical  miles. 


V  =  /OO  /tAfO  r~3 


y  -=  /ooo  y*s. 


Then: 


V  -  r  e 

V  -  Velocity 

r  •  Radius 

h  -  Rate  of  change  of  © 
V 


e  - 


^£    -  200  Radians/iour 


O  =r 


//OVA! 


/  Hoy*  /  £/?£.      _        £  7.  3  4£~<j. 


3£O0  J£c 


/G  b^c. 


/&  3fc . 


&    ~     .  OS&S    **£/^ec 
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Individual  components  within  the  computer  itself  have  accuracies  generally 
of  the  order  of  $   minutes  or  less.   If  a  steady  state  error  to  a  velocitv 
input  of  the  same  mapnitude  is  specified,  then  from  the  relationship: 

This  would  represent  nearly  ideal  tracking  under  the  most  severe  condit- 
ions likely  to  be  incountered  in  a  submarine  fire  control  problem  with- 
in the  next  decade,  and  is  a  value  toward  which  it  would  be  desireable 
to  work.  Obviously,  for  the  great  majority  of  applications  K  would 
not  have  to  be  nearly  as  large  to  obtain  5  minutes  or  less  tracking 
error. 
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V.  Mathenatical  Analysis  of  the  System  and  Comparison  with  Observed  Data 

5.1  Analysis  of  Individual  Loops 

Since  the  four  individual  servo  loops,  Gyro  Anple,  Enabling  Run,  Weapon 

Time  to  Burst,  are  essentially  identical,  only  one  loop  will  be  considered 

in  detail.  This  is  justified  by  measurements  made  on  the  other  loops. 

For  illustration  let  us  consider  the  Gyro  Angle  Loop.  As  previously  stated 

in  section  U.l,  the  pain  cannot  be  considered  in  the  conventional  manner. 

The  closed  loop  equivalent  circuit  is  shown  in  figure  18  below: 
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Fig.  18  CLOSED  LOOP  EQUIVALENT  CTRCUIT  OF  INDIVIDUAL  LOOP 

It  is  noted  that  conventional  summing  resistors  are  not  necessary.  This 
can  best  be  understood  by  considering  the  equivalent  circuit  of  the  input 
and  feedback  vernistat  in  fipure  19. 
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See  Arpendix  I 
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Fig.    19   EQUIVALENT  CIRCUIT  OF  INPUT  AND  FEEDBACK  VERNISTAT 


Further  simplifying: 
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Since  the  reference  voltage  is  variable  and  R       ^   ^        R  the  circuit  can 
be  arranged  in  this  manner. 


-2/r-O.    (Z 


2X0z 


L*:*J 


1 


7~a  /9/#/>jt,/^/&* 


n 


It  is  convenient  to  derive  the  transfer  function  of  the  circuit  comprising 
the  a^lifier  and  integral  damping  network,  figure  20^  prior  to  proceeding 
with  the  analysis. 
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Fig.  20  EQUIVALENT  CIRCUIT  OF  AMPLIFIER  AND  INTEGRAL  DAMPING  NETWORK 

Breaking  the  circuit  in  this  manner  is  justified  since  neglecting  the 
shunting  effect  of  60K  and  R^  -  3  -negohms  results  in  less  error  than 
incurred  in  assumed  values  of  components. 

The  resultinp  block  diarram  is: 


h8 


Simplifying: 
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Fca  ■  closed  loop  transfer  function  of  the  Amplifier. 
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rhe  Root   I/dcus   Eouation  is: 
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(5.2) 


/</       ~    ArooT   *.o<-^j    G^/aj 


^he  close   loop  transfer  function  of  the  amplifier  network  will 


thpr 


have  the  form: 
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The  equivalent  circuit  of  the  entire  loop  then  becomes: 
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The  block  diagram  representation  iei 
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By  block  diagram  manipulation: 
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This  simplifies  to  fipure  21. 
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Fig.  21  REDUCED  BLOCK  DIAGRAM  OF  INDIVIDUAL  LOOP 
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The  root  locus  equation  is: 
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Simplifying: 

5"  (  S+28J*S  t-«J  i5+bJ 

The  closed  loop  function  will  have  the  form: 


£= 


_         20A.  /-#2. 


(S+cj{5  +<tns  +<-  hs+tj 


It  can  be  seen  from  the  above  equations  that  varying  ILj  and  R?  varies 
the  location  of  the  zero  due  to  the  tachometer  and  simultaneously  Ry 
varies  the  root  locus  gain.   It  can  be  seen  that  within  limits  many 
settings  of  R2  and  R-  will  give  identical  results. 

Since  R?  and  Ry  enter  into  gain  as  well  as  zero  location,  the  conventional 
method  of  considering  a  piven  pole-zero  configeration  and  varying  the 
gain  to  determine  optimum  response  cannot  be  used.  Further,  if  a  comp- 
osit  open  loop  transfer  function  is  obtained,  it  can  be  seen  that  these 
factors  also  enter  into  the  computation  of  the  Velocity  Constant,  Kv, 

N 
Usinp  the  relationship  F  «  -** —   the  open  loop  function  is: 

c    c 


Fa   = 


o 


A,= 


(2^t*2)l2 *f*<>  ir^i^trr 


$2 


[t  is  seen  that  three  equations  involving  the  two  qualities,  response  and 
relocity  constant,  can  be  easily  written. 


X     ( Ho*T    Xocuf    GA/AfJ    =      22o*  2ok 


{uok  r/inj 


o      -  22  oo    (jlokj  , 

*7  "  r * 


oK 
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ft-]  =r    ~  A        —       -^  a  A 


(jLotrt-ft) 


«?    =    «*****-  *"   L  W,  S~t°*      +  2j   220  f*   j 

*  -  zvs-x     ~ At*  £    28     +      2  8  fz 

2ok  J 

Since  the  system  cannot  be  optimunrized  by  varying  gain,  a  plot  of  these 
three  equations  will  provide  a  criteria  for  choosing  values  of  resistance  for 
optimum  response  and  maximum  velocity  constant. 

From  figures  22  and  23  it  can  be  seen  that  a  near  optimum  selection  of 
resistance  is  ^*  =  3*°*  and  /f^  -  /***"  .  It  is  noted  that  these 
values  give  Kv  in  excess  of  the  minimum  derived  in  lw3. 
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With  these  values  fixed,  consider  the  loop  in  three  conf igerations : 

1.  Single  loop  -  Unity,  tachometer,  and  integral-damp  feedback 

2.  Sinrle  loop  -  Unity,  and  tachometer  feedback 

3.  Single  loop  -  Unity,  and  integral  feedback 

5#2  Si*Ele  I*QP  -  Unity.  Tachometer,  and  Integ^l-damp  Feedback 

Substituting  R2  -  300  K  ,  R?  .  100  K  ,  and  R6  -  3  Megohms  equation  (5.1) 
becomes: 

/J  =  S5~Q,ooo   iS^ZJ  (9  /-2ooJ 

*  %S  +2j  (5  /-Zoo/    y-    /^r^ys" 


Equation  (5.2),  the  root  locus  equation  of  the  amplifier  and  Integral- 
damp  feedback  network  becomes: 
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The  resulting  block  diagram  is: 
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Equation   (5.3)  becomes: 

(AT,  Joe  J  &+*>  (S  +  2.00J _ 

Simplifying: 

C      s+  +  '?3t5    +  '*2'  ++z  5**  3,Z  *Q,  COOS  /-  t,/zo,ooc 

Equation  (5»U)»  the  final  root  locus  or  characteristic  equation  of  the 
loop  is: 


-/  - 


5(   Si-^8t<Sf-/3J<.J^^.zc^ 


The  root  locus  plot  is  shown  in  figure  2li.  The  roots  are  determined  from 
K  versus  the  real  roots  in  figure  25. 

The  real  roots  were  determined  to  be  S  ■  -  1862  and  S  -  -  2,0332  •  From 
the  characteristic  equation,  denominator  of  equation  (5.b),  the  coeffi- 
cent  of  the  second  highest  term  is  1931  and  the  constant  term  is 
6,120,000.  Proceeding  as  outlined  in  section  3.3: 
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The  factored  closed  loop  function  then  becomes: 
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The  response  to  a  unit  step  input  is: 


Gi.t) , 


/&.3O0    (S/-S7 (&+2O0J 


(S+2.0*3ZXS+'e*&lS  +  3*+j22'&SJ(s+*3-J**-<& 


& 


The  inverse  Laplace  transform  yeilds: 

_  -  Z.0552  *  _   s&izt: 

Ooi*)~    1+       o/a     €  +  .oos    € 

+     /S-2S~      6~S*       S/"(22.6££-/S7.S*J 


Figure  26  shows  predicted  0  as  a  function  of  time.  Figure  27  shows 

the  measured  response  obtained. 

* 

3o2  Single  Loop  -  Unity  and  Tachometer  Feedback 
Removing  the  integral-darnp  network,  the  block  diagram  is: 


9xt) 


J +  .OZ595 
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The  closed  loop  transfer  function  is: 


/:  - 

Simplified  becomes: 


/S.30O 


s(  S  +2$)  /■  J£7  (  s+  */.<*; 


/3~,  30  O 

Sz     +3<fS  S   +     /s:  300 


The  root  locus  eouation  is: 

3£7(S+*'<<> 


-/ 


Si  >t 2,8/ 


The  root  locus  is  shown  in  figure  28.  The  Root  Locus  Gain  versus  Real 
is  shown  in  figure  29. 
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The  factored  closed  loop  transfer  function  is  : 


/=c    = 


/J7  30  o 


{S  +  V¥J    i^/-  3f-sj 

The  response  to  a  unit  step  is: 

/S,3J°  (4-) 

The  inverse  Laplace  transform  yields s 

-  fV  £ 
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The  predicted  output  as  a  function  of  time  is  shown  in  figure  30.  The 
measured  response  is  shown  in  figure  31« 


3.3  Single  Loop  -  Unity  and  Integral-damp  Feedback 
With  no  tachometer  feedback,  the  block  diagram  becomes: 
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The  block  diagram  simplifies  to: 
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The  closed  loop  transfer  function  then  becomes: 


K» 


SiT,  3QO    (  S+ZJ  ** +2J0' 


Si  S+ZMS+.21J  (5  T-/S-JL;   /■ 


/s;  *o°  <  s+Z/i? *-^3<*j 


6$ 


66 


s 


H 

W 


o 

M 

En 


o 

a 


£3 


w 
o 


o 
o 


a. 
w 

Eh 
CO 


W 

CO 

is 

o 
a. 


rt 


CO 


M 

=3 

E- 


67 


Simplified: 


£*  s*+,sm*'  +*».**'***   7  3.,~.»o.s   *  *.'*:"<> 


The  root  locus  equation  is: 


/s.304  (5/- 2 j  (sy-^ao/ 


=  -1 


The  root  locus  plot  is  shown  in  figure  32.  The  plot  of  K  versus  Real 

Roots  is  shown  in  figure  33 •  The  factored  closed  loop  transfer  function  i3: 


r   _  /S,joo     <JS+1)  Ury-z-oo/ 


(St^  Z- 6+1/(5  J~  /S2-V/& +  /<•>«■*  +/  <*i*8A  S+/i'Hfy  *■/•#/ 

The  response  to  a  unit  step  is  piven  by: 


The  inverse  Laplace  transform  vields: 

r/  y.  .  a/iVs  c      +  ■  oas  7S  e 
**"[*<■  /.oz  £-*•"' S,~  iW.9*-  **.fj  J 

The  predicted  output  as  a  function  of  time  is  shown  in  figure  3U.  The 
measured  response  is  shown  in  figure  35. 

5.2  Qualitative  Stability  and  Performance  Analysis  of  Overall  System 
Section  3.1  outlines  the  Determinantal  Method  of  analysis  used  in  reducing 
the  functional  diagram,  figure  36,  to  the  configeration  shown  in  figure 
37c 

It  is  noted  that  the  Torpedo  Lead  angle  Loop  is  inherently  non-linear 
because  of  the  resolver  which  is  used  for  a  function  generator  and  a 
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System  Functional  Diagram 


nulling  point.  The  output  from  the  resolver  is  Zl  Cos  Bb6  -  XX  Sin  Bb6 
which  tends  to  zero  when  Bb6  is  computed  correctly. 

Because  of  this  non-linearity,  it  was  assumed  that  the  inputs  Z   X  Cos  Bb6 
and  Z  T  Sin  Bb6  were  constant  inputs  and  not  a  function  of  time.  The 
resolver  was  represented  as  a  summer.  Several  of  the  actual  constants 
were  combined  with  the  scale  factors  of  the  lxjop   for  convenience  and  simp- 
licity. This  is  justified  since  it  does  not  effect  the  accuracy  of  the 
stability  analysis  and  it  is  highly  unlikely  that  quantitative  response 
can  be  predicted  with  the  assumptions  made  concerning  the  resolver  unless 
a  prohibitive  number  of  small  perturbations  are  considered. 

The  assumption  is  also  made  and  justified  in  section  6.1  that  for  sipnals 
which  do  not  saturate  the  amplifier  or  motor  that  the  Enabling  Run, 
Weapon  Time  to  Burst,  and  Gyro  Lead  Angle  transfer  functions  are  essentia- 
lly identical.  For  large  signals  this  assumption  is  especially  invalid 
in  the  Weapon  Time  to  Burst  Loop  because  of  the  I4O  to  1  gear  ratio  which 

permits  the  motor  to  saturate. 
1 

With  these  assumptions  in  mind,  the  "Standard  Block  Diagram"  becomes 
that  shown  in  figure  37.  Where: 
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Looking  at  the  output: 
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It  23  convenient  to  expand  the  delta  so  that  the  J\    ~  can  be  obtained 
since  it  will  be  used  in  the  numerator  of  the  closed  loop  function. 
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Simplified 

A   =  /  r    G,G.G7    r  ftft  ft  *  G  ft  ft  ft-ft  *  #  &  c-  ft 

The  closed  loop  function  becomes 


/-     ^  &-  /,   6.  if  4.     63    J  J      *    Wf    U*.    <*«    J 


/    /    (7,  <?„  £7    *  CI  ft  C*  /■  <£-£  G,  C*  G,  Cs  +  G.  G<  G¥j 


Using  the  relationship  /^ 


ASc 


0*-Mc. 


the  open  loop  function  is 


/- ,  ft/lftft  gB  ,*  ftflftj 
/  ^  ft  ft  ft  *  <?„  ft  <?<, 

Since  each  function  G.  C2.      ®tc.  has  the  form  of  1   numerator  divided  by  a 


denominator,  let  us  define 


Gj    =    j£L 


Substituting  and  simplifying: 


L~  „        M  "±  Afj  ASs-  A/q     7*      AS,  Ar^  A/*  AS*-    0j    £><.    &?  &t 

&  />*  £>*  0<  £>j-  A  #t/>9   ?  /^  "*  *7  &  2>j  £>4-4e/>e  *>i  At,  *L  A  4  ^-  ^  #s 

In  terms  of  numerators  and  denominators  of  the  various  functions,  the  closed 

loop  function  becomes: 

f  =     ^ *  1  At,  My  Afj  ^v  ^  0^     +    Af.  A*x    /*V  A/f  0j  0t  &T0* 

0.^0,  Jv  *r4c  #7  09+  -*i  "v"7  *<  ^  -4- A^    "**i  A"t  4  0+  24-*r  4 
A   >K  y^*  Afj  AS^A-g  0*  A  ^7    *     ^  >*t  *V  >K-  ^  ^c  A  <^ 
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Substituting  the  various  values  for  the  numerator  and  denominator  of  the 
functions,  the  open  and  closed  loop  transfer  functions  become: 

US+IUD*  {5+  2.0332J*  (S+33+f  22.LS)^\ 
[(JS  +  33 -J  22.  CS^  J 

f  (5  *  \V.l)*-(  5>+  2.0S3Z)*-  (Si-    S3  -/•/  22  4*7* 
\j-   Ur.300)*  (5+2J*LSi'*0eJ*    /fx   (  *-+/&+*,+,] 


-         VS139*)*  *.k*\.  kik+-*s-k(.\(S+xJ*(J>+2ooJ¥ 

[(   .S>  /7(.2)K  (Si-  2.  0331)*-  (S  +  33  +j  22.  C*J  \| 
"(  $r  /6L2)X  I  S  +2.0  332)*  (  SV-  33  +j  22.  C  3') 
+  Ux  3*o/*"  (5**;*  C3  +2.co;*-   fx  (*+*j-  +  trg  kj 

+  {/S;  3oo)  *t  ks  l**kH~  *>  kj    LS+2J+(S+Jco)+ 


Observing  equation  5.5  the  open  loop  transfer  function  is  of  the  form 

F  •  *<*+*'' l***i J    where  n  -  0.  This  indicates,  of  course,  that  the 
system  is  type  zero. 


(SSJ 


-  ^tUw*  s 


[k'  (S+£.j  is**»; 1 
(S+4H&+/U J 


The  steady  state  error  to  a  ramp  is: 

3 


£  * 


This  indicates  that  the  steady  state  error  would  be  infinite  and  the  sy- 
stem must  be  compensated  in  such  a  way  to  make  it  a  type  one  if  ramp 
inputs  of  large  magnitude  or  long  duration  are  expected  to  be  incountered. 
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Since,  of  necessity,  certain  linerizations  were  made;  therefore,  no 
conclusions  can  sensibly  be  made  concerninr  stability  except  that  the 
necessary  condition  that  all  the  coefficients  of  the  characteristic 
equation  are  positive  is  met.  This  being  the  circumstance,  a  frequency 
response  analysis  was  made.  The  results  were  plotted  on  a  Nichol's 

chart.   It  was  found  that  there  was  a  gain  margin  of  26  db  and  a  phase 

o 
margin  of  72  .  This,  of  course,  shows  that  the  overall  system  is  highly 

stable.  The  results  obtained  from  the  Nichols  chart  were  then  plotted 

on  a  Bodie  diapram,  figure  38 »  to  determine  the  dominant  time  constant 

of  the  system  and  to  acertain  experimentally  that  the  system  was  actually 

type  zero.  These  results  are  discussed  in  section  5.3-b. 

5.3  Discussion  of  Analysis  and  Observed  Data 

5»3a  Individual  Loop 

The  individual  loop  was  analyzed  in  three  possible  configerations, unity, 
tachometer,  and  integral-damp  feedback;  unity  and  tachometer  feedback; 
and  unity  and  integral-damp  feedback.  * 

By  using  figure  22  and  23  it  was  found  that  near  optimum  values  of 
resistance  were  Rp  ■  300  K  and  R?  ■  100  K.  It  was  found  that  the  value 
of  R/  used  had  little  effect  on  the  response  of  the  system;  however, 
considerations  of  noise  in  the  system  dictated  a  value  of  approximately 
three  megohms.  Using  these  values  of  resistance,  computations  were  made 
to  predict  the  response  of  the  loop.  The  predicted  and  measured  response 
is  shown  in  firure  26  and  27  respectively. 

Prior  to  the  determination  of  the  values  of  resistances,  the  best  response 
obtained  was  a  critically  damped  case  with  a  settling  time  in  excess  of 
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,2h   seconds.  It  can  be  seen  that  the  predicted  response  and  measured 
response  agree  exactly  and  a  settling  time  of  .1$  seconds  was  obtained. 
This  shows  an  improvement  in  response  of  37.5  per  cent. 

The  analysis  of  the  loop  without  the  integral-damp  network  was  performed 
in  section  3.2  •  The  measured  response  is  shown  in  figure  31.  Comparing 
this  with  the  predicted  response  in  figure  30,  it  is  again  seen  that  the 
measured  and  predicted  values  agree.   It  is  also  noted  that  an  improvement 
in  response  was  obtained  with  a  settling  time  of  .11  seconds.  This, of 
course,  is  desireable;  however,  at  steady  state  there  was  audible  noise 
caused  by  gear  chatter.  This  noise  did  not  show  up  in  the  output  but 
would  cause  excessive  wear  in  the  gear  train  near  the  motor;  therefore, 
the  integral-damp  network  is  considered  desirable.  It  was  found  analy- 
tically that  by  changing  capacitors  C^  and  Ci  (See  figures  13  or  lii)  to 
,2$M-f  a  predicted  and  measured  response  with  a  settling  time  of  .12 
seconds  was  obtained  without  any  associated  noise  problems.  The  theo- 
retical pole  location  is  calculated  in  section  b.2. 

The  last  configeration  considered  involved  unity  and  integral-damp  feed- 
back. From  figures  3U  and  35  it  can  be  seen  that  the  predicted  and 
measured  response  are  again  in  agreement  except  the  peak  overshoot  is 
somewhat  greater  in  the  actual  system.  Settling  time  is  .3  seconds. 
It  is  noted  that  substantial  tachometer  feedback  must  be  used  since  for 
smooth  computation  each  loop  must  be  critically  damped. 

5.3b  Overall  System 

As  discussed  in  section  5.2,  the  overall  system  is  highly  stable  with  a 

o 

phase  margin  of  72  and  a  gain  margin  of  2k   db.  The  data  from  the  Nichol's 

plot  is  used  to  construct  figure  38.  From  this  it  is  noted  that  the 
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first  asvmptote  has  a  zero  db  slope;  therefore,  type  zero  system.  The 
time  constant  agrees  reasonably  with  the  actual  obtained  from  the  meas- 
ured response  shown  in  figure  39.  The  theoretical  settling  time  is  .6 
seconds  and  the  measured  i3  .7  seconds  (See  figure  39).  The  response 

to  a  ramp  input  is  shown  in  figure  UO. 

» 
Since  the  theoretical  computations  and  freauency  response  indicated  that 

the  system  was  type  zero  with  a  finite  error  to  a  step  input  and  an 
infinite  error  to  a  ramp  input,  it  was  considered  advantageous  to  meas- 
ure these  errors.  The  inputs  and  resulting  errors  are  shown  in  figures 
Ul  and  U2. 

The  results  are  therefore  conclusive  that  the  system  is  in  fact  type  zero. 
While  the  multi-loop  system  is  unmodified  (i#e.,  a  type  0  servomechanism), 
there  will  be  a  constant  error  to  a  step  and  a  corresponding,  increasing 
error  to  a  ramp  input.  The  measurements  in  figures  Ul  and  U2  show  the 
magnitude  of  the  errors  involved.  The  first  is  a  measurement  of  the  error 
versus  time  for  the  step  and  the  second  is  the  measurement  of  the  error 
versus  time  corresponding  to  the  ramp.  From  the  first  the  error  ratio 
to  the  step  input  can  be  calculated  to  be: 

( 
Tn  a  typical  submarine  problem,  the  total  bearing  change  from  initial 

contact  to  the  firing  point  is  approximately  1$0  degrees.  This  means 

that,  under  these  circumstances,  the  error  at  firing  point  will  be: 

One  degree  «  17.6  yards/1000  yards.  Therefore,  the  following  errors 
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would  result  with  a  computed  error  in  the  Gyro  Lead  Angle  Computation 


Loop  of  C.8  : 


att  1,000  yds. 
2,000  yds. 
3,000  yds, 
li, 000  yds. 


Hi. 05  yds. 
28.10  vds. 

U2.15  yds. 
56.20  yds. 


These  errors  are,  of  course,  in  addition  to  any  other  errors  generated 
in  other  parts  of  the  commuter.   It  can  be  readily  seen  that  if  fire 
control  problems  of  short  range  with  the  resultant  large  ramp  or  long 
ranges  where  tracking  is  done  for  long  periods  of  time,  a  significant  error 
may  d eve lope. 

Looking  at  equation  ($. 5>), the  open  loop  transfer  function  of  the  system, 
it  is  seen  that  an  integration  is  required  or  an  "Sw  is  needed  in  the 
denominator  to  have  a  type  one  system.   It  appears  that  a  feed  forward 
loop  added  to  the  block  diagram,  figure  37,  would  accomplish  the  desired 
results.  The  resultant  standard  block  diagram  is  shown  in  figure  U3. 
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Fig.  U3  STANDARD  BLOCK  DIAGRAM  WITH  COMPENSATION 
NECESSARY  TO  PRODUCE  A  TYPE  ONE  SYSTEM 
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The  determinant  for  delta  becomes 

a         fc>        c 
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o     *f     O      ' 


o      o 

O       O 

O        O      -Cs      -G+     ,        q 
Gf      0       O        O       o     -C/-     /    I 

The  determinant  reduces  to 

A  =    /  /   <?*  <?v  <?7    +  G>.GsC*    +  <?f    +  <?*  <?v  C7  Of    *  <7,  Gs  0.  G9 
+    G<  <?*  <?v  Gs   +  <?,  G.  G,  Cr  Cs 

Aa+  -     C7  /  Gz  (7V  07  G?  +  Gz  Gs  <?<  CV  *  tf  <T*  <7V  C^  /■  (^  <?,  £,  C<J> 

The  closed  loop  function  is 

[/*<■  Gz  <T*  G7  /  Gz  Gj  Gc   *  Gf  /■  £  (?«,  cT7  c^   *<?*  <?,  C«  CV 
1/  C  <?*  CV  Cr   *  G.  Gz  G,  (?r  G, 

In  terms  of  numerator  and  denominators 


J 


+    A/<  Af*  Af?  AsrAQ  0*  0c  4h  &f 

A  Pz  A/>y  A-  te  #7  P*  #r     +  ^  "+  *i   *'  &   *r  ^  ^7  ***  4r 
/    Mi.  **  Afc    />,/>«  />r  J7&9  #7  +  S*9AJ>z  0jP«  P^AcPyft 

+   Ss±  AV   *7  *<;   f<  />J    4r  ft  ^   +  /Vl  #*  A*c  **r  4*  **  ^^^7  ^ 

/    jr,/lsi    /Is*  "4-  />j  Z>c    #7  &6  &r    *  /U'  **  ^^  ><V^V64^ 
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The  open  loop  function  is 

/**   /).  Vx  6>i  Js*  &s-  &<.  Vi  *>8    +   /ft-  #*  *7  *i  &  &  £r-#<.    &9 

■/■     A/±  #3  yf«.  As*    #'    f+   £j-  #1  #8     /•   A/***.  /*>  #r  *>*  4*.  £7   £*  "f 

IT  +  AS*  ^    /*- j   /*>  /^y  ^>*  />*  />?  /V 

/>/  />*  ^j  />y  ^J»  /4  ^7  ^  /^   *  /^*.  /*V  X>  ^  0j#r-6><.   4  ^  ^ 
/■  /</*  /V,  /^«,  />,  /V  At  //7  ^  ^> 

It  can  be  seen  that  the  necessary  and  sufficient  condition  for  a  type  one 
system  has  been  met  since  fi9  »   3C  S+&J      can  be  factored  from  each  term. 
Where  G.  -  *1   is  of  the  form:  G^  -  ^ff^J 
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71.  Conclusions  and  Recommendations 

6.1  Single  Loop 

The  single  loop  designed  and  constructed  by  Librascope  Division  of  General 
Precision  Corporation  will  produce  satisfactory  response  without  further 
changes  in  design.  However,  by  judicious  selection  of  the  various  input 
resir.tors,  an  improvement  is  possible  of  37.5  percent  using  settling 
time  as  a  criteria.  It  can  be  seen  from  the  root  locus,  figure  2\\   and 
K  versus  Real  roots,  figure  2£  that  the  system  is  absolutely  stable  for 
all  values  of  gain.  It  can  also  be  seen  that  the  gain  can  vary  appreci- 
ably without  effecting  the  response.  Should  the  gain  become  very  low, 
an  appreciable  residue  will  be  associated  with  the  root  near  the  zero 
at  S  ■  -2.  Tachometer  feedback  is  essential  to  provide  a  critically 
damped  system  for  computation.  Integral-damp  feedback  slows  down  the 
response  by  26.5  per  cent  from  that  obtained  without  it;  however,  the 
damping  is  required  to  prevent  high  frequency  oscillations  in  the  gear 
train  during  steady  state  operation.  In  selecting  values  of  capacitance 
(Co  -  C»  -  .25><"r  )  in  figure  Ik,   the  determental  effect  on  response 
due  to  the  integral-damp  feedback  is  eliminated.  The  recommended  changes 
improves  the  response  a  total  of  5U.2  per  cent  over  the  best  obtained 
by  trial  and  error  on  the  physical  system. 

6.2  Overall  System 

The  overall  system  is  type  zero  as  designed.  This  is  satisfactory  for 
most  applications  intended  since  the  error  involved  is  small;  however, 
should  large  ramp  inputs  or  tracking  problems  of  long  duration  be  encount- 
ered, it  is  possible  to  generate  excessively  large  errors.  The  solution 
to  this  problem  is  a  modification  involving  a  motor  and  amplifier  in  a 
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feed  forward  loot).  This  loop  will  be  physically  connected  at  the  error 
input  to  the  Torpedo  Lead  Angle  and  feed  forward  to  the  output  of  Bb6. 
This  is  represented  in  the  modified  standard  block  diagram,  figure  h1. 

An  alternate  method  which  may  be  prefered,  since  the  primary  function 
of  the  loop  is  to  mechanize  an  equation,  is  to  determine  a  mean  ramp, 
calculate  the  various  outputs  as  a  function  of  time,  and  scale  the  vari  us 
loops  to  reduce  or  eliminate  the  resultant  error,  A  check  must,  of  course, 
be  made  over  the  full  range  of  ramps  expected,  after  the  scaling  has  been 
accomplished,  to  insure  that  the  error  does  not  become  excessive  else- 
where. 

Measured  frequency  and  transient  response  indicates  that  the  system  operation 
will  be  satisfactory  with  .7  seconds  settling  time  to  a  unit  step,  gain 
margin  of  2h  db,  and  phase  margin  of  12   . 
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VII.  Recommendations  for  Further  Investigation 

Should  it  be  deemed  desirable  to  predict  accurately  the  response  to  a 
complex  multi-loop  system  utilizing  resolvers,  it  would  be  advantageous 
to  use  the  linearizations  made  herein,  consider  small  pertubations,  and 
determine  their  accuracies  obtained  bv  measurements  on  the  physical  system. 
This  would  provide  a  criteria  for  evaluating  future  analysis  which  may  be 
conducted  during  design  work  on  other  systems. 

In  addition  it  may  be  advantageous  to  consider  linearization  technioues 
other  than  the  ones  used  herein. 

If  it  is  reauired  that  the  overall  system  be  optimized,  an  analogue  co"nputer 
study  is  warranted. 
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APPENDIX  I 

Induction  Resolver 

An  induction  resolver  is  an  electromagnetic  device  employed  as 
a  basic  element  in  analog  computers,  automatic  control  systems, 
angle  data  transmission  systems  and  plan  position  indicating 
radars.  It  is  used  for  performing  trigonometric  computations. 

The  induction  resolver  is  essentially  a  transformer  with  rotary 
variable  coupling  between  primary  and  secondary  windings  (Figure  1). 
The  windings  are  distributed  in  the  various  stator  and  rotor  slots 
in  such  manner  as  to  obtain  a  true  sinusoidal  relationship  between 
output  voltage  and  angle  of  rotation  of  the  rotor.  Since  the 
units  are  used  at  audio  and  power  frequencies,  it  is  customary 
to  use  a  high-permeability  laminated  core  similar  to  that  used  in 
an  induction  motor. 


^=/;  s/^y 


Fig.  1  -  Basic  Resolver 
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It  is  customary  to  treat  the  stator  winding  as  the  primary  or 
input  winding,  and  the  rotor  winding  as  the  secondary  or  output 
winding.  If,  as  in  Figure  2,  there  are  two  secondary  windings 
at  right  angles  to  each  other  (in  "space  quadrature"),  one  sec- 
ondary will  have  induced  in  it  a  voltage  which  is  proportional 
to  the  sine  of  the  angle  of  rotation  between  it  and  the  primary, 
and  the  other  secondary  will  have  induced  in  it  a  voltage  which 
is  proportional  to  the  cosine  of  the  angle  of  rotation  between 


it  and  the  primary. 
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Fig.  2  -  Resolver  with  two  Rotor  Windings 
To  make  induction  res^lvers  even  more  useful,  two  windings  are 
provided  on  both  the  stator  and  rotor.  The  pairs  of  stator  wind- 
ings and  rotor  windings  are  located  at  right  angles  to  each 
other  as  shown  in  Figure  3.  This  permits  computation  of  cosines 
as  well  as  sines,  and  even  more  complex  problems  involving  both 
sines  and  cosines  can  be  solved  when  both  stators  are  excited. 
In  addition  to  the  resolution  of  input  voltages  into  sine  and 
cosine  components,  resolvers  can  perform  the  synthesis  of  resul- 
tants from  components  and  the  rotation  of  coordinate  systems. 
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Fig.  3  -  Re solver  Input  and  Output  Relationship 

As  shown  in  Figure  h,   currents  in  the  two  stator  windings 
set  up  magnetic  fields  at  right --angles  to  each  otrer  to  form  a 
resultant  flux  vector.  This  flux  vector  may  assume  any  angular 
position,  depending  upon  the  amplitudes  and  polarities  of  the 
input  currents.  Its  magnitude  and  proportional  to  the  vector 
sum  of  the  instantaneous  values  of  the  input  currents,  and  its 
angular  position  is  equal  to  the  arctangent  of  the  ratio  of 
the  two  currents.  For  AC  input  voltages  within  the  operating 
frequency  range  of  the  resolver,  the  rotor  voltages  are  propor- 
'-.tional  to  the  magnitude  of  the  flux  vector  and  the  sine  or  cosine 
of  the  angle  of  rotation  of  the  rotor  winding  with  respect  to 
the  vector.  The  output  voltages  may  be  expressed  in  terms  of 
the  input  voltages  and  the  rotor  angular  position,  by  the  fol- 
lowing equations: 
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f„^e5     Co*Q-£i     J>/sVG  (1) 

The  subscripts  refer  to  the  various  rotor  and  stator  windings 
and  &  is  the  rotor  angular  position. 


...     ,   ! ^  '  «-«- 

Fig.  h   -  Stator  Flux  Vectors, 

While  other  types  of  operations  are  possible  with  a  resolver, 
its  use  has  been  principally  with  three  types  of  inputs,  namely; 
(1)  sine  wave  voltages  of  the  same  electrical  phase,  (2)  sine 
wave  voltages  shifted  90  degrees  apart,  and  (3)  irregular  wave 
shapes,  principally  of  the  sawtooth  or  square  wave  variety.  For 
the  first  type  of  operation,  the  output  voltages  remain  fixed 
in  electrical  phase  shift  and  their  magnitudes  vary  in  accordance 
with  equations  (1)  and  (2).  For  the  second  and  third  types  of 
operation,  the  output  voltages  may  also  be  determined  by  equations 
(1)  and  (2)  if  the  instantaneous  values  of  the  voltages  are  taken 
into  consideration.  For  the  case  of  the  resolver  with  two  inputs 
90  degrees  apart  in  electrical  phase,  each  of  the  two  rotor  output 
voltages  remains  constant  in  amplitude,  independent  of  rotor  an- 
gular position.  However  the  phase  of  the  voltages  varies  contin- 
uously with  rotor  angular  position.  The  phase  shift  between  the 
two  output  voltages  remains  fixed  at  90  degrees.  In  this  mode  of 
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operation,  then,  the  equations  for  the  output  voltages  are: 

^ .=    e,   IS.  (3> 

For  the  third  case,  where  irregular  wave  shapes  are  used,  op- 
eration most  frequently  falls  into  two  categories.  In  phase, 
square  waves  or  flat  topped  sine  waves  a^e  used  in  some  computing 
applications.  The  performance  is  adequately  described  by  equations 
(1)  and  (2)  without  modification  if  the  two  input  wave  shapes 
are  the  same,  which  is  usually  the  case. 

Sawtooth  wave  shapes  are  frequently  used  in  radar  sweep  aypli- 
cations.  In  the  simplest  case,  a  sweep  voltage  is  applied  to 
one  stator  and  the  two  output  voltages  from  the  rotor  windings 
are  the  same  sweep  voltage  but  with  an  amplitude  proportional 
to  the  sine  or  cosine  of  the  rotor  angular  position.  If  these 
two  voltages  are  applied  to  the  nX"  and  "Y"  deflection  plates, 
or  coils  of  a  cathode  ray  tube,  they  cause  a  radial  trace  to  be 
formed.  The  trace  will  assume  the  same  angular  position  as  the 
resolver  rotor.  This  type  of  operation  finds  application  in 
plan  position  indicator  types  of  radar  displays. 

While  the  arrangement  shown  in  Figure  3>a  is  suitable  for  schema- 
tic representation  of  a  resolver,  tie  arrangement  shown  in  Figure 
5b  is  more  suited  for  block  diagrams  and  functional  schematics. 
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The  diagram,  equivalent  schematic  and  equations  are  shovn  for 

clockwise  rotation  (facing  the  shaft  end)  and  apply  to  all 

standard  types  of  resolvers.  Vhere  it  is  desired  to  change  the 

direction  of  rotation  of  the  resolver,  it  is  only  necessary  to 

interchange  input  A  with  B  and  output  C  with  D.  This  can  be 

noted  -^n  the  block  diagram  by  a  negative  sign  preceding  q   . 
A 
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A  resolver  used  in  the  simplest  form  of  computation  is  shown  in 
Figures  6a  and  b.     One   stator  winding  is  excited  with  a  voltage 
proportional  to  the  hypotenuse  of  a  right  triangle  and  the  rotor 
is  positioned  to  the  angle  <^    .     The  two  output  voltages  are  pro- 
portional to  the  two  legs  of  the  triangle. 
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Fig.  6  -  Elementary  Resolver  Application 
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In  Figure  7  the  problem  is  similar,  but  the  hypotenuse  and  one 
of  the  legs,  Y  are  known.  A  servo  amplifier  and  motor  are  used 
to  position  the  resolver  rotor  until  the  voltage,  ^R«,  is  del- 
ivered by  the  r««olver  equals  the  known  voltage  Y.  The  voltage 

v 
R_  is  then  proportional  to  the  length  of  the  other  leg  of  the 

triangle,  and  the  angle  ©  is  the  angle  between  the  second  leg 
and  the  hypotenuse.  A  solution  of  this  type  may  be  used  in  de- 
termining elevation  angle  and  ground  range  when  slant  range  and 


depth  are  known. 
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Fig.    7  -  Angle  Computation 

In  Figure  8,  the  inputs  are  two  vectors  90  degrees  apart.     The 
servo  positions  the  resolver  rotor  until  one  output  is  zero. 
At  this  position  the  other  rotor  being  90  degrees  away,  has  maxi- 
mum coupling  with  the  flux  vector  formed  by  the  X  and  Y  inputs. 
Its  output  is  therefore  proportional  to  the  vector  sura  of  these 
two  inputs.     The  rotor  angular  position  is  the  angle  between  one 
leg  and  the  hypotenuse  of  the  triangle. 
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Fig.   8  -  Vector  Addition 
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In  the  rotation  of  coordinates,  both  stators  and  rotors  are 
■used.  V.lth  rectangular  coordinates  X  and  Y  as  inputs  to  a  resol- 
ver,  as  in  Figure  9,  the  outputr-  are  also  rectangular  coordinates 
of  the  same  point  referenced  to  a  new  set  of  axes  rotated  with 
respect  to  the  first.  The  operation  if  fully  described  by  the 
equations  in  Figure  3  which  will  be  recognized  as  those  required 
for  coordinate  rotation. 


The  actual  operation  can  best  be  understood  by  breaking  it  into 
two  steps.  The  X-Y  coordinates  impressed  on  the  two  stators 
form  a  flux  vector  in  the  re solver.  This  vector  is  independent 
of  rotor  position.  Its  angular  position  is  the  arctangent  of  Y/x. 
Each  of  the  rotor  windings  develops  a  voltage  proportional  to 
the  sine  or  cosine  of  its  angular  position  with  respect  to  the 
flux  vector  (not  with  respect  to  zero).  If  the  reference  is  made 
with  respect  to  zero  then  the  outputs  follow  the  equations 
in  Figure  9.  ^e  rotor  voltages  are  the  X',  Y'  components  of 
the  flux  vector.  The  angle  between  the  X*,  Y'  coordinate  ref- 
erence of  the  inputs  i§  the  same  as  the  rotor  angular  position. 
V. 
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Fig.  9  -  Coordinate  Rotation 
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INDUCTION  POTENTIOMETER 

An  induction  potentiometer,  or  linear  synchro  transmitter,  may- 
be considered  a  special  kind  of  resolver  which  provides 
accurate  linear  indication  of  shaft  rotation  about  a  reference 
position  in  the  form  of  a  polarized  voltage  whose  magnitude  is 
proportional  to  angular  displacement,  and  whose  phase  relation- 
ship indicates  direction  of  shaft  rotation.  The  principal  differ- 
ence between  a  resolver  and  an  induction  potentiometer  is  that 
the  latter' s  output  voltage  varies  directly  as  an  angle  and  not  as 
the  sine  of  that  angle.  Its  rotation  is  usually  limited  to  less 
than  +90°. 
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Fig.  1  -  Induction  Potentiometer  Voltage  Output 
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These  devices  are  analogous  to  resistance  potentiometers,  but 
since  they  are  induction  type  components  they  have  less  re- 
straining forces  acting  upon  their  rotors,  and  hence  are  capable 
of  providing  better  resolution.  Naturally,  induction  poten- 
tiometers do  not  require  sliders  such  as  those  used  in  resistance 
types.  Therefore,  circuit  interruptions  are  eliminated,  no  wear 
occurs  as  a  result  of  rubbing  parts,  and  accuracy  is  consequently 
continously  maintained  at  the  original  level  throughout  the 
operational  life  of  these  types  of  components.  These  advan- 
tageous features  are  virtually  a  necessity  in  certain  applications. 
For  example,  in  gyroscope  systems  where  low  restraining  torques 
and  low  pick-off  angular  errors  are  required,  wide  use  is  made 
of  these  devices.  In  general,  induction  potentiometers  find  use 
in  applications  where  resistive  potentiometers  are  impractical, 
principally  because  of  the  following  features:   (1)  induction 
potentiometers,  having  no  wiping  contacts,  may  be  used  as 
gyroscope  pick-offs  since  they  contribute  less  spurious  friction 
torque:   (2)  input  and  output  are  isolated:   (3)  resolution  is 
infinite:   (U)  noise  level  is  low:  and  (3>)  the  total  angle  of  travel 
is  limited  to  less  than  180 

Constructed  very  much  like  a  resolver,  induction  potentio- 
meters differ  in  that  their  windings  and  slots  are  not  uniformly 
distributed,  but  instead  are  deliberately  modified  to  produce  a 
linear  output.  Normally  they  have  only  one  excitation  winding 
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and  one  output  winding,  the  former  usually  carried  on  the  rotor 
and  connected  by  means  of  two  sliprings  and  two  brushes. 

The  use  of  induction  potentiometers  as  linear  pick -of fs  may 
be  considered  as  a  special  resolver  application  in  which  mutual 
coupling  varies  not  as  ^S/ir^,  but  as  f  9«  With  aero  input  im- 
pedance  and  infinite  output  impedance,  the  voltage  gradient, 


fer  function  may  be  expressed  as: 


■  .  Under  this  condition  the  trai 
as: 


H,  may  be  defined  as  H  •  — .  Under  this  condition  the  trans- 


where  &»    •  secondary  impedance 
*>  »  primary  impedance 
j?a  «  load  impedance 
Z-ji  "  input  impedance 
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